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in ' ABSTRACT 

(N ' 

q ■ We study a combined sample of 264 star-forming, 51 composite, and 73 active 

■ galaxies using optical spectra from SDSS and mid-infrared (mid-IR) spectra from 

■ the Spitzer Infrared Spectrograph. We examine optical and mid-IR spectroscopic 
diagnostics that probe the amount of star formation and relative energetic con- 
tributions from star formation and an active galactic nucleus (AGN). Overall we 

5a ■ find good agreement between optical and mid-IR diagnostics. Misclassifications 

of galaxies based on the SDSS spectra are rare despite the presence of dust ob- 
scuration. The luminosity of the [Nell] 12.8 jttm emission-line is well correlated 
with the star formation rate (SFR) measured from the SDSS spectra, and this 
holds for the star forming, composite, and AGN-dominated systems. AGN show 
a clear excess of [Nelll] 15.6 \xm emission relative to star forming and composite 
systems. We find good qualitative agreement between various parameters that 
probe the relative contributions of the AGN and star formation, including: the 
mid-IR spectral slope, the ratio of the [NeV] 14.3 fim to [Nell] fim 12.8 fluxes, 
the equivalent widths of the 7.7, 11.3, and 17 \xm PAH features, and the optical 
"D" parameter which measures the distance a source lies from the locus of star 
forming galaxies in the optical BPT emission-line diagnostic diagram. We also 
consider the behavior of the three individual PAH features by examining how 
their flux ratios depend upon the degree of AGN- dominance. We find that the 
PAH 11.3 fim feature is significantly suppressed in the most AGN-dominated 
systems. 



*S. M. LaMassa is now at Yale University. 
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Introduction 



It has been established that a strong link exists between supermassive black holes 



(SMBHs) and the galaxies in which they live (IFerrarese fc Merrittl |2000| ; iGebhardt et al. 



2000l ; iTremaine et al.l 120021 ) and betwee n accreting SMBHs (or active galactic nuclei, AGN) 



and star formation in particular (e.g. iKauffmann et al.l 120031 ; ICid Fernandes et al.l 1200 ll ) 



Clues to this link can be revealed by detailed analysis of the spectra of galaxies and AGN. 
Star formation processes imprint their signatures on the spectrum. Stellar populations of 
various ages produce strong continuum emission in the ultraviolet (UV) through near-IR 
region, while gas photo-ionized by the hot young stars produces prominent emission-lines 
in the UV, optical, and IR. The accretion disk surrounding the supermassive black hole in 
an AGN prod uces its own strong UV and optical continuum emission (the "big blue bump 
( IShieldslll978f )). This emission can swamp that of the stellar population of the host galaxy 
when there is a direct view of the accretion disk (i.e. in Type 1 AGN). Dust in the interstellar 
medium of the host galaxy will absorb optical and UV light from both stars and the AGN 
and re -radiate this in t he IR. In AGN, this re processing also occurs in the dusty obscuring 
torus (|Antonuccilll993t lUrry &: Padovanilll995l ). In type 2 AGN, the direct view of the AGN 
is blocked by the obscuring torus, and the presence of the AGN can be inferred by prominent 
optical and infrared emission lines from highly ionized gas located beyond the torus and the 
strong mid-IR continuum emission from the torus. In such objects the signatures of host 
galaxys young stars and star formation are also present in optical and mid-IR spectra. These 
type 2 AGN are therefore ideal laboratories for studying the connection between AGN and 
star formation, provided that the tracers of this activity can be disentangled. 

Multi-wavelength star formation indicators have been extensively studied in samples 
of quiescent star forming galaxies (see lKennicuttlll998l . for a review). Star formation rates 
(SFRs) have been calibrated in the optical based on Ha emission, which results from recom- 
bination following ionization due to energetic photons from massive O and B stars, and in the 
UV based on the continuum which provides a window into emission from young stars. Dust 
in the galaxy reprocesses optical and UV photons and re-radiates in the infrared (IR), thus 
revealing obscured star formation. Other IR indicators include emission from fi ne structure 
lines excited by starburst activity, such as [Nell] 12.81/im and [Nelll] 15.56/im fjHo fc Keto 
20071 ). as well as spectral emission features associated with polycyclic aromatic hydrocar- 
bons (PAHs), since these macro - molecules /small- grains can be ionized by single optical or 
UV photons dLi fc Draind 120021 ) and are prominent in IR spectra of star-forming galaxies 
f lSmith et alll2007bh . 



Such host galaxy star formation diagnostics can be contaminated by the presence of an 
AGN. Photons emitted by the accretion disk enhance the optical and UV emission used to 
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derive SFRs. The obscuring medium enshrouding the central engine contributes significantly 
to the infrared emission, co nstituting; approximately 20% of the AGN bolometric luminosity 
( jSpinoglio fc Malkanlll989f ). Disentangling the contribution from the AGN vs. star formation 
then becomes a necessary intermediate step in using continuum IR emission as a tracer of star 
formation in active galaxies, and such a process can be plagued with uncertainties. Previous 
studies have indicated that emission from PAHs tend to be suppressed in galaxies hosting 
an active nucleus, which is sometimes attribu ted to the harsh radiation field produced by 
the AGN accretion disk destroying: PAHs (e.g. |Voit|ll992l; b'Dowd et alJ l2009l: llrever et ah 



2010l ; lLaMassa et "aDboiol : iDiamond-Stanic fc Riekelboioh . 



A complementary view of the AGN and starburst connection is to study spectroscopic 
signatures that parameterize the relative importance of these two processes. The optical 
BPT diagram (a plot of [OIII]A5007/H/3 vs. [NII]A6583/Ha) provides a useful diagnostic 
to differentiate between star-forming galaxies, composite systems (galaxies with comparable 
amounts of star-formation and AGN activity) and Type 2, or obscured, AGN in the local 



universe ([Baldwin et al.l Il98ll ). The optical "D" parameter is the distance a source lies 



from the locus of star for ming galaxies on the B PT diagram: a higher D value indicates 



greater AGN dominance (IKauffmann et al.l 120031 1 . Ratios of IR fine structure lines also 



primarily ionized by A GN (IRigby et al.ll2009l ; IDiamond-Stanic et al.ll2009t lArmus et al 
Goulding et al. 2010 ) whereas [NeII]12.81/im is excited by star formation ( Ho fc Ketol 2007 ). 
The ratio of these l ines can then indica t e the relative importance of these two processes 
(IGenzel et al.l Il998l ; lLaMassa et al.l |2010| ; iPereira-Santaella et al.l |2010| ) . The mid-infrared 
(MIR) spectral slope, 020-30^, steepens as the amount of emission from cold dust heated 
by stars increases relative to the amount of emission from hotter dust heated by the AGN. It 
is thus anoth er potentially useful tool to assess the rel ative amount of star formation in active 
systems (e.g. iBuchanan et al.l 120061 ; iDeo et al.ll2009l ). The equivalent width (EW) of PAH 
grains is another diagnostic to probe the relative amount of star formation to AGN activity. 
An anti-correlation between ionization field hardness and the ratio of PAH features at 7.7/xm 
and 11.3/xm has been reported and is interpreted as due to an increasing contribution to the 
mid-IR continuum emission by AGN-heated dust. The selective destruction of smaller PAH 
grains may also play a role (e.g. ISmith et alll2007bl ; lO'Dowd et al.ll2009t iTreyer et al.l 12010 



Wu et alJEoiol ). and so the flux ratios of the different PAH features can then also potentially 
trace the relative contribution of AGN and star formation in active galaxies. 

In this paper, we combine samples of star-forming galaxies, obscured AGN (Seyfert 
2 galaxies, or Sy2s) and composite systems to study the interplay between AGN activity 
and star formation. Using quiescent star forming galaxies as a base-line, we investigate the 
effects of AGN activity on the following star formation diagnostics: optically derived SFRs 
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from the Sloan Digital Sky Survey (jBrinchmann et a 



structure lines [Nell] 12.81/xm and [Nelll] 15.56/um fjHo fc Ketol 120071 ) and the luminosities 



2004), the luminosities of the IR fine 



of the polycyclic aromatic hydrocarbons (PAHs) at 7.7/im , 11.3/xm and 17/im fjSmith et al 



2007bl ; iFarrah et al.l 120071 ; iDiamond-Stanic fc Riekd l2010l ). Which of these star formation 



proxies agree the best among star-forming galaxies and AGN and are thus least affected by 
the presence of an AGN? 



We expand upon the work of lLaMassa et al.l (120101 ) to test the agreement among di- 
agnostics that parameterize the relative contributions of AGN activity and star formation, 
i ncluding the equivalent width (EW) of PAHs at 7.7, 11.3 and 17/im, ratios of PAH fluxe s 



JSmith et all l2007bl; loT towd et al.lbood: iTrever et al.lboiol: biamond-Stanic fc Riekelboioh 



ratio s of IR fine structure lines (iGenzel et al 



20101 ). mid-IR sp ectral slope (a^o- 



1998; 



Armus et a 



Buchanan et al 



cal D parameter ( Kauffmann et al. 20031 ). In 



200G; 



LaMassa et al 



112007; Pereira-Santaella et al. 



Deo et al.ll2009i ) and the opti- 



(boioh we found that Sy2s with 



stronger PAH emission tend to have a softer ionization field and that the MIR spectral slope 
was well correlated with PAH EW. By expanding this parameter space into the regime of 
quiescent star forming galaxies, we test if PAH features and the mid-IR spectral slope are 
dependent on the hardness of the radiation field, regardless of its source, or if ionization by 
an AGN is necessary to appreciably affected observed features. Finally, we use the results of 
this analysis to empirically decompose the mid-IR (MIR) emission into a star-forming and 
AGN component. 



2. Sample Selection 



To assess the efficacy of star formation diagnostics in active galaxies, we have com- 
bined samples of AGN, quiescent star for ming galaxies, and c omposite objects. Sources 
were categorized by their BPT designations (iBaldwin et al.lll98ll ) using the boundaries from 



Kewlev et all (hoOlf ) and iKauffmann et al.l (bo03h : galaxies above the iKewlev et al.l (120011) 



delineation are c 



Kauffmann et al 



Kauffmann et al 



a ssified as pure Seyfert 2s, sources in between the IKewlev et al.l (120011 ) and 
( 2003 ) demarcations are classified as composites and the sources below the 
(120031 ) boundary are classified as quiescent star-forming galaxies. 



As discussed in lLaMassa et al.l (120091 ). we selected an optical sample of Sy2s from the 
Sloan Digital Sky Survey (SDSS). These sources were identified as pure AGN on the basis 
of their [OIII]A5007/H/3 and [NH| A6583/Ha flux ratio s which place them above the theo- 
retical starburst demarcation from lKewley et al.l (120011 ) on the BPT diagram. This sample, 
complete to an [OIII] flux limit of 10~ 14 erg/s/cm -2 (hereafter "[OIII] sample"), comprises 
20 Sy2s. We complement these Sy2s with a mid-infrared (MIR) selected sample, the orig- 
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inal 12/xm sample ISpinoglio Sz Malkanl ( 119891 ) from the IRAS point source catalog, which 
contains 31 galaxies. Using published or SDSS measured values of [OIII]A5007/H/3 and 
[NII]A6583/Ha, we have classified 27 12/xm sources as Sy2s, 2 as composites (F08572+3915 
and NGC 7590) and 2 were dropped due to lack of published [Nil] /Ha values (F04385-0828 
and NGC 3982). As both [QUI] and MIR fluxes are indicators of intri nsic AGN power (e.g 
Spinoglio fc Malkanl Il989t iHeckman et al.l 120051 ; lLaMassa et al.l 120101 ) , these Sy2 samples 
should to first order represent unbiased views of the AGN population. The [OIII] sample 
has optical spectroscopy from SDSS whereas the optical data for the majority of the 12/im 
sample were based on litera ture values. Both samples have low and high resolution Spitzer 
Infrared Spectrograph (IRS iHouck et al.l 120041 ) data. 



Our comparison samples of predominantly normal galaxies have SDSS spectroscopy, 
Spitzer IRS data and ultraviolet photometry from the Galaxy Evolution Explorer (Galex). 
The Spitzer-SDSS- GALEX Spectroscopic Survey (SSGSS, P.I. Schimi novich) contains 101 



local (0.03 < z < .21) galaxies oyerlap ping the Lockman Hole region (jO'Dowd et al.ll2009 



Treyer et al.l l2010t lO'Dowd et al.ll201ll ). This sample is augmented by the Spitzer SDSS 
Spectroscopic Survey (S5, P.I. Schiminovich), increasing the sample size of local (in this 
case, 0.05 < z < 0.1) star-forming galaxies to 300 (Bertincourt et al., in prep). Both the 
SSGSS and S5 samples have IRS Short-Low (SL) and Long-Low (LL) coverage. The S5 
sample also has IRS Short-High (SH) for all sources, but only 1/3 of the SSGSS sources (a 
sub-sample consisting of the brightest galaxies) have high resolution infrared spectra. Based 
on the BPT designations, there are 26 Sy2s, 49 composites and 264 star-forming galaxies in 
the SSGSS and S5 surveys. 

Combining all four samples, we garner 264 star-forming galaxies, 51 composites, and 73 
AGN for this analysis. Throughout this paper, we have adopted a cosmology of H Q = 70 
km/s/Mpc, n M = 0.27, and Q x = 0.73. 



3. Data Analysis 
3.1. Optical Data 

Optical spectra for the [OIII], SSGSS and S5 samples were drawn from SDSS, using the 
MPA/JHU galaxy value-added ca talog@ This catalog includes ancillary derived data, such 
as star formation rates (SFRs, iBrinchmann et al.l 120041 ). which we utilize in the subsequent 
analysis. Only four of the 12/im selected sources have SDSS spectra (3 Sy2s and 1 composite). 



2 http:/ /www. mpa-garching.mpg.dc/SDSS/DR7/ 



-6 - 



Optical data fo r the r emainder of these sources were culled from the literature, as noted in 



LaMassa et al.l (120 lOh . Hence, these 24 Sy2s and 1 composite are omitted from the SDSS 



SFR analysis in Section 4. We emphasize here that all SDSS-derived quantities are based 
on the SDSS spectra alone. 



3.2. Infrared Data Reduction 



in 



The details o f the S pitzer data reduction f or the [QUI] and 12/im sam ples were presented 



LaMassa et al.l ( j2010h . the SSGSS sample in lO'Dowd et al.l (120091. 120111) and the S5 sampl e 



in Bertincourt et al., in prep., which follows the same procedures as lO'Dowd et al.l ( 120091 ). 
In short, the IRS data were run through the standard Spitzer calibration software. Sky 
subtraction was performed for all low-resolution data and for the high resolution sources that 
had dedicated off source observations. All but 68 observations, from the SSGSS sample, have 
high resolution data, leaving 223 star-forming galaxies, 33 composites and 64 Sy2s with high 
resolution IRS coverage. Only eight sources, from the 12/x sample, did not have dedicated 
high resolution background observations, but since we utilize these data to measure emission 
line fluxes, the addition of the background to th e underlying continuum has negligible effect 
on our derived results (see lLaMassa et al.ll2010l ). S pectra were extract ed with SPICE for the 
[OIII], SSGSS and S5 sources and with SMART rtffigdon et alJ liooih for the staring mode 
12/im observations. CUBISM fjSmith et al. Il2007al ) was used to extract spectra from 12/im 
spectral mapping observations. 



3.3. Emission Line Measurements 

Emission line fluxes used in this analysis ([Nell] 12.81/mi, [NeV] 14.32/im and [Nelll] 
15.56/ahi) were measured from the high resolution spectra, specifically the Short-High (SH) 
module. The local continuum around the emission feature was fit with a local zero or first 
order polynomial. This local continuum fit was frozen and the emission feature was then fit 
with a Gaussian. The error was estimated by smoothing the local continuum by the line 
width (or the nominal minimum resolvable line width§ in the case of non-detections) and 
calculating the RMS around this smoothed continuum. In order for a line to be considered 
detected, it had to be significant at the 5a level or above, have a line width at least on the 
order of the SH resolution (i.e. full width half max > 0.01 fim) and have a line centroid 



3 width = A nom /R, where X nom is the nominal wavelength of the emission line and R is the resolution of 
the SH module, ~600 
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near the nominal wavelength of the emission line (i.e. < 0.03/im). The [Nell] and [Nelll] 
fluxes and upper limits for the SSGSS, S5, [OIII] and 12/xm samples are listed in Tables [1] 
- HI Since we consider [NeV] upper limits in aggregate rather than individually (see Section 
7.1), we list only the [NeV] detections in these tables. 



3.4. PAH Measurements 



We measured the fluxes and equivalent widths (EWs) of the polycyclic aromatic hy- 
drocarbons (PAH) at 7.7/im, 11.3/xm and 17/im with PAHFIT, a spectral decomposition 
routine that models dust continuum, starlight con tinuum, unresolved emission lines, PAH 
emission features and silicate absorption/emission ( ISmith et alll2007bl ). PAHFIT requires a 
single stitched low-resolution spectrum. The Long-Low (LL) spectra was therefore shifted 
to match the Short- Low (SL) spectra in their overlapping regio ns, with typical adjustments 
~20% or lower (see lO'Dowd et al.ll201ll ; lLaMassa et al.ll2010l ). This choice of normaliza- 
tion was made since the SL and SH apertures sizes are well-matched, better enabling us to 
compare quantities derived from the high and low resolution spectra. 

The routine was run with a mixed extinction model on all sources. We note that in 



LaMassa et al.l ( 120101 ) . PAHFIT was run with a screen extinction, which assumes a fore- 
ground absorbing slab rather than dust mixed with gas, starlight and grains. We have 
re-run PAHFIT on the [OIII] and 12/im sources using the default mixed extinction option 
to be consistent with the SSGSS and S5 samples, though differences in results were slight. 
For these two Sy2 sa mples, we also a dded a dust component at 1000K since this may better 
model strong AGN ( iDeo et al.l 120091 ). though we note the difference in the fluxes and EWs 
between the default temperature array and this modified version is negligible for a majority 
of the sources and a factor of 2 or lower for most of the remaining sources. We consider 
a PAH feature to be detected if it is present at greater than the 5<r level. In the subse- 
quent analysis, we consider PAH upper limits on the flux, but only include EWs for detected 
sources. PAH measurements are listed in Tables - El 



4. Parameterizing Star Formation: Effects of AGN Activity 



Optical and IR fluxes have been utilized to calibrate SFRs in quiescent galaxies. Below 
we investigate if emission from the AGN contaminates these parameters and if composite 
systems are more similar to star-forming galaxies or Seyfert 2s. Here we cons i der op tical 
SFRs derived from SDSS using the methodology described in iBrinchmann et al.l ( 120041 ) . We 
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also explore IR star formation proxies, namely the fine structure [Nell] 12.81/xm and [Nelll] 
15.56/im emission lines as well as the luminosity of the polycyclic aromatic hydrocarbon 
(PAH) features at 7.7, 11.3 and 17 /xm. 



4.1. SDSS SFR 



fiber • 



Circumnuclear Star Formation 



For star-forming galaxies, SDSS SFR/j(, er is calculated from the Ha luminosity, corrected 
for dust extinction, measured through the 3" SDSS fiber. In these systems, Ha emission 
results from the recombination of ionizing photons from massive O stars , making this flux a 
usefu l parameterization of recent star formation. Us ing a Kroupa IMF (IKrqupa fc Weidner 
20031 ). Ha luminosity is then converted into a SFR (IBrinchmann et al.ll2004l ). 



Such a straightforward calculation fails for active systems due to emission powered by 
the AGN contaminating the Ha flux. However, a SFR can be estimated in these systems 
using the strength of the Dn(4000) feature, which is a break in the optical spectrum at 4000A 
that indicates the age of the host galaxy stellar population. The SFR/jb er for composites and 
AGN in SDSS is calculated starting with the strength of the Dn(4000) break. This is then 
converted into a SFR by adopting the relationship between SFR_f,;h Rr /M» , nhe.r and Dn(4000) 



meas ured for quiescent systems and then using measured value of M*jj& er (IBrinchmann et al. 
20041 ). 



4.2. IR SFR diagnostics: Neon Emission Lines and PAH Luminosities 



In the infrared, [Nell] at 12.81/mi has an ionization potential of 21.56 eV and is therefore 
excited by star formation activity. [Nelll] at 15.56/im has a higher ionization potential than 
[Nell] at 41 eV, b ut it become s more dominant in lower-metallicity star forming systems (e.g. 



Wu et all 120061 ). iHo fc Ketol (120071 ) have demonstrated that the sum of [Nell] and [Nelll] 



correlates better with total infrared luminosity (a reliable proxy of SFR) than just [Nell] as 
this sum better accommodates the low metallicity systems. They subsequently calibrated a 
SFR based on [Nell] + [Nelll]. 

PAH macro-molecu les /small- grains are ionized by single ultraviolet or optical pho- 
tons (ILi fc Draind 2002l ) and are very prominent in the spectra of star forming galaxies 



( Smith et all 



Farrah et al. 



2007bf) . SFRs have consequently been calibrated on PAH luminosities (e.g. 



20071 ). However, the PAH luminosity may be affected by AGN act i vity, and the 



effect may be different for the diferen t PAH features (e.g. lO'Dowd et al.ll2009t iTreyer et al. 



20101 ; iDiamond-Stanic &: Riekei |2010| ). In this analysis, we consider the sum of the PAH 
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luminosities at 7.7, 11.3 and 17/xm (Lspah7. 7,n.3,i7^m)- 111 Section 5, we examine each PAH 
feature individually as a function of increasing ionization field strength. 



4.3. Comparison of Optical and IR SFR diagnostics 



We compare the optical and IR SFR diagnostics pair- wise in Figures [TJ- [3] as both scatter 
plots (left hand panels) and histograms of the ratios of the two proxies under consideration 
(right hand panels) to search for systematic offsets among galaxy typeQ In these plots, 
blue represents star-forming galaxies, green denotes composites and red indicates AGN. 
To take into account upper lim i ts, we used survi val analysis (ASURV, iFeigelson fc Nelson 



19851 ; llsobe and Feigelsonl Il990l ; llsobe et all Il986l ) to determine the best-fit trends to the 



star-forming galaxies and Sy2s (shown as the blue dotted-dashed lines and red dashed lines, 
respectively, in the left hand panels of Figures [T|- [3j) , to calculate the mean ratio of SFR prox- 
ies for each galaxy type (Table [9]) and to quantify the significance of agreement /disagreement 
of the distributions of the SFR proxy ratios among star-forming galaxies, composites and 
AGN (Table PTOl) . We use the Gehans Generalized Wilcoxon Test to calculate the probabil- 
ity that the distributions of SFR proxy ratios among the sub-samples being considered are 
drawn from the same parent population. A probability value (P) < 0.05 indicates that we 
can reject the null hypothesis and that the distributions differ at a statistically significant 
level. We note that not all sources are included in each calculation, either due to absence of 
data (i.e. 25 12/xm sources have no SDSS coverage) or lack of significant detections in the 
diagnostics being considered. The number of sources used for each calculation are listed in 
Tables M -ECU 

As shown in Figured] and Table ITU} the [Nell] luminosity in relation to SDSS SFR^ er . 
is consistent among star-forming galaxies, composites and AGN. However, the sum of [Nell] 
and [Nelll] does show signs of AGN contamination as these sources are systematically offset 
to higher L[ NeII ^ + L^ NeII n values with respect to star-forming galaxies, as seen in Figure 
[21 Indeed, the distribution of Lr Ne jn + L^em] normalized by SFRjj 6er is significantly dif- 
ferent between star-forming galaxies and AGN (Table [TO]) . This result is due to the harder 
ionization field from the AGN boosting [Nelll] to higher values, biasing the sum as a SF in 



dicator in active systems. This result is consi s tent with previous works (e.g.lDeo et al. 



Tommasin et al. 2008: Melendez et al. 2008; Wu et al. 



20m IPereira-Santaella et all 120101 ) 



2007 



4 In these plots and subsequent plots of PAH fluxes and EWs, the error bars are smaller than or on the 
order of the symbol size and therefore not plotted. These relatively small errors are due to the stringent 5tr 
threshold we required for detection. 
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that demonstrated that [NeIII]/[NeII] is correlated with AGN activity and thus a tracer of 
ionization field hardness. Interestingly, the composite galaxies are not affected: their Ljjve/j] 
+ Lrjvem] sum differs significantly from AGN, but not from star-forming galaxies. Hence, 
the SFR calibrated on Ljjve/j] and ^[Neiii] is reliable for star-forming and composite galaxies, 
but not AGN. Consequently, to test the efficacy of PAH luminosities in parameterizing star 
formation in AGN, we compare PAH luminosities to SDSS SFR/j;, er and L[jv e /r| an d not the 
sum of L[jvei7] and L [NeIII] . 

Both [Nell] and PAH luminosities are upper limits for at least a handful of sources. 
Since ASURV is unable to accommodate doubly censored data points, we calculated the ratio 
means and performed two sample tests using upper limits on [Nell] and PAH luminosities 
separately. We report the results of both analyses in Tables M and [TD1 including the number 
of sources for each calculation. We find that the distribution of PAH luminosities relative 
to SFR/jfe er . is consistent among star- forming galaxies, composites and AGN, though the 
slope of the relation for Sy2s differs markedly compared with star-forming galaxies (Figure 
[3]). The difference in PAH emission becomes more pronounced when compared to Lrjv e /n, 
regardless of whether L^p^y 711 317Mm or Ljjve/j] is treated as the upper limit (see Table [TU]) . 
As Table M demonstrates, there is a decrease of about 32% in the mean in the ratio of PAH 
to [Nell] luminosity in AGN compared to both the quiescent star-forming galaxies and the 
composites. At face value we can not say whether this is caused by supression of the PAH 
emission or enhancement in the [Nell] emission. We will return to this issue later, when we 
explore the effects of ionization field hardness on the relative strengths of PAH features at 
7.7, 11.3 and 17 pm individually. 



5. Relative Importance of AGN to Star Formation Activity 



Expanding on the work of lLaMassa et all (120101 ). we now investigate the agreement 
among multi-wavelength diagnostics of ionization field hardness to determine where "nor- 
mal" star-forming galaxies live in this parameter space. As AGN activity weakens and star 
formation becomes more dominant, do the trends found for AGN smoothly extend into the 
regime of quiescent galaxies? 



As discussed in lLaMassa et al.l (120 lOf ). the [O IV]26/xm/[Ne 



able tracer of the incident radiation field hardness (Genzel et al. 



2010 ) since [OIV] is primari l y ionized by AGN acti vity ( Rigby et alibooijl : 



2009 



I]12.81/xm flux ratio is a reli- 



19981: 



Pereira-Santaella et al. 



Diamond-Stanic et al. 



Melendez et al. 



2008 



LaMassa et al.ll2010l ) whereas [NeII]12.81/xm is excited by star 



formation (IHo fc Ketol 120071 ). However, we do not have Long-High (LH) coverage for the 



SSGSS and S5 samples, and thus lack high resolution data for the [OIV] line at 26/im. 
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We therefore use [NeV] /[Nell] as a diagnostic of the relative imp ortance of AGN versus 
starburst activity ( lArmus et al.ll2007t iPereira-Santaella et al.ll2010l ). We compare this ratio 



with other diagnostics studied in lLaMassa et al.l (120 lOh . First we use the optical D param- 
eter (D = y/[log{[NlT]/Ha) + 0.45] 2 + [log([OIII}/H(3) + 0.5] 2 ), which is the distance a 
source lies on the BPT diagram from the locus of star- forming galaxies, so tha t a higher D 
value indicates a greater dominance of AGN activity (IKauffmann et al.ll2003l ). Second we 
use a 20 -3o M m (MIR spectral i ndex) which increases i n the presence o f cold dust associated 
with star formation activity ([Buchanan et al.l 120061 ; iDeo et al.l 120091 ). Finally we use the 
PAH equivalent widt hs (EWs) which a re expected to decrease as AGN activity becomes 
more dominant (e.g. iGenzel et al.lll998l ). We also investigate if the ratios of the fluxes o f 



the different PA H fe atures are affected by AGN activ ity, as found by lO'Dowd et al.l (120091 ). 



Wu et all (l2010l ) and lDiamond-Stanic fc Riekel (J2010|). 



5.1. [NeV] 



Due to its hi gh ionization potential (97.1 eV), [NeV] has been claimed to clean ly trace 



SMBH accretion (IGenzel et al.lll998l ; lArmus et al 



be an unambiguous signature of AGN activity (jGilli et al. 



2007 



Pereira-Santaella et al.l 20101) and to 



2010 



Goulding et al.ll2010l ). Also 



as it is formed in the narrow line region, it should ideally be unaffected by toroidal obscura- 
tion a nd therefore a tracer of intrinsic AGN flux, s i milar to the [OHI] 5007A and [OIV] 26/xm 



2009 



LaMassa et al 



lines ( Heckman et al.l 120051 ; iMelendez et al.l 120081 ; iRigby et al.l 120091 ; iDiamond-Stanic et al. 



2010|). However, in a majority of our spectra, this feature was not 



significant above the noise. To use [NeV] /[Nell] as a tracer of ionization field hardness, we 
consider the non-detections in aggregate, stacking the spectra among sub-groups of galaxy 
type, rather than employing individual upper limits. For this stacking exercise, we separated 
star-forming galaxies, composites and AGN into two groups each, based on their location in 
the BPT diagram, in order of increasing nuclear activity. Figure H] illustrates the method we 
used. The l eft hand plot is a BP T diagram for our full sample, with th e unbroken black lines 
showing the iKewley et al.l (1200 ll ) (upper) and IKauffmann et al.l (120031 ) (lower) demarcations 
between star-forming galaxies, composites and AGN. The dashed lines illustrate the addi- 
tional empirical separations we used, approximately tracing the halfway point within each 
region. These lines correspond, in order of lowest to highest, to the following equations: 



Log([OIII\/Hp) 
Log([OIII]/Hf3) 



0.61 



Log([NII]/Ha) + 0.1 
0.61 

Log([NII}/Ha) -0.2 



+ 1.3 



1.3 
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LoMOIIiyHB) = Logm ™ l /Ha) _ , + 1.19 

We also employed a cut on [Nil] /Ha to exclude low metallicity systems: Log ([Nil] /Ha) 
> -0.6 dex for star-forming galaxies and Log ([NII]/Ha) > -0.2 dex for AGN. The right- 
hand panel of Figur e H] shows ou r resulting sub-groups d epicted in different colors with the 



Kewley et al.l (120011 ) (upper) and iKauffmann et all (120031 ) (lower) demarcations overp lotted 



star-forming group 1 - dark blue, star-forming group 2 - light blue, composite group 1 - dark 
green, composite group 2 - light green, AGN group 1 - gold, AGN group 2 - red. 

We stacked together the IRS SH spectra for the sources in each color bin from Figure H] 
b) not detected in [NeV] . Each spectrum was shifted to the rest-frame, rebinned to a common 
wavelength grid of equal pixel spacing and normalized by the the average MIR flux between 
13.25/im and 13.75/xm; this window was chosen since it is free of strong emission features. 
Figure \5\ shows the stacked SH spectra for each galaxy bin. We used a 3a clipping on the 
average flux in each pixel to reject outliers that would have otherwise introduced erroneous 
features into our stacked spectra due to noise. [NeV] was detected above the 5a level in just 
one of these bins, AGN group 2. For the other bins, we obtained [NeV] upper limits using the 
procedure outlined in Section 3. [NeV] /[Nell] was calculated by normalizing the [NeV] upper 
limit /detection by the stacked [Nell] emission feature in each group, which was detected at 
high significance in all cases. In the following analysis, we plot the [NeV] /[Nell] detections 
and stacked [NeV] /[Nell] upper limits against the optical D value, PAH EW (for the sub-set 
of sources where the PAH luminosity was detected above the 5a level) and 0:20-30/^ • For the 
stacked [NeV] /[Nell] ratios, we calculate the average D, PAH EW and a20 — 30/xm value in 
each bin. The [NeV] stacked detection is illustrated by the triangle in the subsequent figures. 

Figures [6] - [7] illustrate the trends between [NeV] /[Nell] and the other diagnostics of 
radiation field hardness. As expected, [NeV] /[Nell] is significantly correlated with D and 
significantly anti-correlated with a2o-30/im and PAH EWs, as seen in the correlation coef- 
ficients (p) and associated probabilities that the two quantities are uncorrelated (P un corr) 
reported in Table [Hj We note that one galaxy classified as a quiescent star-forming galaxy 
based on its optical (SDSS) spectrum has a [NeV] detection. This suggests that it is har- 
boring a "hidden" AGN, i.e. an AGN embedded in very dusty host galaxy that attenuates 
the optical signatures used to identify AGN. However, we note that this source has other 
IR parameters (i.e. 020-30^™, PAH EWs) more consistent with those of typical star-forming 
galaxies. This may argue against "hidden" AGN interpretation, which may then call into 
question using solely [NeV] detection as an unambiguous signature of AGN activity. 

In any case, we emphasize that the optical and mid-IR classifications of galaxies as 
being dominated by star formation or by a Type 2 AGN agree in the overwhelming majority 
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of cases. This is good news, since it shows that SDSS classifications (which are available 
for roughly one million galaxies) are highly reliable. We do not find a large population of 
optically hidden AGN, as found by Goulding & Alexander (2009). Over half of the [NeV] 
detected sources in their volume limited sample of nearby galaxies are not optically classified 
as AGN. However, they do not include optically classifed LINERS or composites as "optically 
classified AGN," while we do consider these to contain AGN. The fraction of [NeV] detected 
sources in their sample that are classified optically as quiescent star-forming galaxies is 24% 
(4 of 17), which is still a substantially larger fraction than what we find. 



5.2. 



Q!20-30/im 



The m id-IR spectral slope index, a?2n-30um, in creases with the relative amount of star 
formation ([Buchanan et al.l 120061 ; iDeo et al.l 120091 ). We plot a2o-30^m as a function of the 
optical D parameter in Figure [H] which shows that a wide range of a 2 o-3o^ m values are present 
for the star-forming galaxies. Consequently, there is no significant anti-correlation present 
in the full sample, with correlation coefficient (p) of -0.195. However, if we consider just 
the AGN, a significant trend is apparent: p=-0.474. This correspondence likely drives the 
significant anti-correlation found for [NeV] /[Nell] vs 020-30^™ • The IR spectral index is thus 
useful in parameterizing the relative contribution of star formation in active galaxies but 
this appears to break down in the quiescent galaxy regime. Consequently, we do not utilize 
c^2o-30fim as a tracer of ionization field hardness in the subsequent analysis. 



5.3. PAH EW and Flux Ratios 

Previous studies have indicated that the PAH equivalent widths (EWs) are anti-correlated 
with the presence of increasing AGN activity. This can be due to two effects. First, as the 
AGN becomes increasingly dominant energetically, the emission due to warm dust heated 
by the AGN will rise. Even if the intrinsic PAH emission is unaffected, this increase in the 
mid-IR continuum flux will cause the PAH EWs to decrease. A second possible effect would 
be due to the actual dest ruction of dust grains by the hard ionization field from the AGN 



central engine ( 1Voitlll992l ). Our new expanded data set that spans the range from quiescent 
systems to strong AGN provides the opportunity to test these processes over a full range of 
ionization field hardness and to test which features may be most sensitive to AGN activity. 

Figure M plots the PAH EW at 7.7 /xm, 11.3 pm and 17/im as a function of the optical 
D parameter. The relationship between D and the PAH EWs is significant, with p <-0.65 
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(Table [T2|) . However, if we focus on just the star- forming galaxies and composites, this 
significant anti-correlation disappears: p=-0.227, -0.125 and -0.230 for the 7.7pm, 11.3pm 
and 17 pm features, respectively. Hence, it appears that only in the objects that are optically 
classified as AGN is the PAH EW appreciably affected, suggesting the AGN contributes a 
significant fraction of the mid-IR continuum or perhaps destroys PAHs. 

As noted above, the PAH EW are sensitive to both the relative amount of mid-IR 
continuum from the AGN as well as effects on the PAH emission itself. A complementary 
approach that isolates the effect of the AGN directly on the PAH emission is to consider 
changes in the ratios of the different PAH featur es as the AGN becomes more domin ant. 
Thus, similar to the work of O'Dowd et al. ( 2009 ) and Diamond-Stanic fc Riekel ( 201ol ). we 
investigate the relationship between radiation field hardness (parameterized by [NeV]/[NeII] 
and D) and the ratios of the PAH fluxes in Figures [TU] and [TTJ Since the [NeV] line is only 
detected in the optically-classifed AGN (with few exceptions), we compare D to the PAH 
ratios for the AGN subset alone, as well as for the full sample (including star-forming and 
composite systems). Qualitatively, it appears that the distribution of the AGN are different 
in these figures from the star forming and composite objects. For the most AGN-dominated 
sources (D > 1.2 and log(L[j\r e v]/L[AreW]) > — 0.5), there is a subset of the population that 
has abnormally low ratios of PAH 11.3 pm to 17 pm and abormally high ratios of PAH 
7.7 pm to 11.3 pm. A similar but weaker effect is seen in the PAH 7.7 pm to 17 pm 
ratios. These distributions do not appear to take the form of a gradual correlation across 
the whole galaxy population, but rather appear as a threshold effect at high values of D and 

l°g (L[7V e y]/L[7v e //]). 

Despite this caveat, we have quantified the statistical significances of overall trends that 
may be present. We emphasize that the correlations may be weak, but the trends found can 
be quantified as statistically significant as indicated by the correlation coefficients (p) listed in 
Table IT3"1 L P AH7 \7 fim /^pahu >m displays no significant trend with ionization hardness, regard- 
less of the radiation field parameterization used and whether the full sample or AGN only are 
considered. Using [NeV] /[Nell] as a probe of the ionization field, a significant trend becomes 
apparent with Lp AH11 , 3tim /Lp AH17fim (p=-0.564) and L PAH7 , 7fim /L PAH11 , 3fim (p=0.611). How- 
ever, with D as a diagnostic of the incident radiation, the Lp AHU3fim /L PAH i 7fim relationship 
becomes less significant in the full sample (p=-0.364) and even less correlated for the AGN 
sub-sampple (p=-0.247). l^pAm.7fim/^-'PAHii.?,fim is not significantly correlated in the full sam- 
ple (p=0.179), but becomes more signficant in the AGN sample (p=0.431). Taken together, 
these results suggest that the PAH feature at 11.3pm is suppressed in AGN-do minated sys- 
tems re l ative to the 7.7 and 1 7 pm features. Th is dis agrees with the results of ISmith et al 
J2007bh . lO'Dowd et all J2009L IWu et all feoioh and biamond-Stanic fc Riekel feoioh who 
find that ratio of Lp A h 7 -. 7 fim/Lp A Hii.3fim is smaller in AGN than in star-forming galaxies. 
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The disagreement with the first three studies listed above could result from the inclusion 
of a greater number of high-luminosity AGN-dominated in our analysis. Clearly, AGN affect 
PAH macro-molecu les and this effect could be more p ronounced in more luminous AGN. The 



inconsistency with iDiamond-Stanic Sz Riekd (120 10[ ) could be due to their Seyfert sample 



being systematically closer. Perhaps the AGN effect on PAHs a re different in the region 



hundr eds of parsecs around the central engine, as probed by the IDiamond-Stanic fc Rieke 



(120101 ) analysis, rather than the kiloparsec to galaxy-wide scale that we are investigating 
where there is a greater interplay between host galaxy star formation and AGN activity. 
Different screening effects at different physical scales may play an important role. More 
advanced theoretical work on the effect of the AGN ionization field on PAH macro-molecules 
would help to understand the physical causes of these inconsistencies. 

Because we have combined a number of different samples selected in different ways, our 
AGN systems (mostly drawn from the SDSS [OIII] and 12/xm samples) are typically at lower 
redshift than most of the star-forming and composite systems (drawn from SSGSS and S5). 
Could any of the results above be influenced by the systematic differences in the physical 
sizes being probed (i.e. smaller scale nuclear emission in closeby AGN versus larger scale 
host galaxy emission for sources further away)? To test this, we impose a redshift cut of 
z>0.03, which corresponds to the lowest redshift source in the SSGSS sample. This cut 
eliminates most of the 12/im sources and almost half of the [OIII] sample. However, since 
we do retain a handful of [OIII]-selected AGN which are strong AGN, we are able to test the 
relationship between PAH EWs and ionization field hardness over a wide range of nuclear 
galactic activity. Figure [12] and [HI] illustrate that the trends found for the full sample still 
hold for this pared down sample. The AGN at higher D values are associated with suppressed 
PAH EWs at a statistically significant level even when the host galaxy is more fully sampled 
(see Table [T2|) . The only change in the PAH ratio analysis is the significant anti-correlation 
of L>pAHii.3iim/L>pAHi7iim with respect to D in the AGN sub-sample (see Table [T3"j) . We note 
that though the L P AH7.7 f im/^-'PAHii.3^m trend is driven by a handful of AGN, the relationship 
is similar to the parent dataset where the parameter space of high LpAH7.7/j,m/^PAHii.3tim 
values is more fully sampled. The thresholds in D and [NeV]/[NeII] for changes in the PAH 
ratios noted above are still present. 

Taken together, these results suggest that some care be used when using the PAH 
luminosities to estimate star formation rates in galaxies with dominant AGN. The 11.3 fim 
feature appears to be rather strongly suppressed in the most AGN-dominated systems. 
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6. Disentangling Starburst and AGN Contributions to MIR Emission 

Mid-infrared emission in active galaxies resul ts from both the toroidal o bscuring: medium, 



accounting for ~20% of the bolometric AGN flux (ISpinoglio fc Malkanlll989l ). and host galaxy 
star formation. We define Lmir as the rest-frame luminosity at 13.5/im averaged over a 3/im 
window: this region is free from strong PAH features and provides a reliable estimate of the 
MIR continuum emission. Using quiescent galaxies as a baseline to estimate the amount of 
MIR emission due to star formation, we can empircally separate the contributions of star- 
burst activity and toroidal emission. To quantify the relative importance of these two effects, 
we solve: 



LmIR — OiLsFR + 0L 



AGN 



We set Lsfr = Lpveir] and Lagn = ^[oiii],agn, the AGN co ntribution t o the o bserved 



[OIII] emission. To calculate L[ ni] ,agn, we use Equation 3 in I Wild et al.l (120101 ). which 
scales the observed [OIII] flux based on the source's position in the BPT diagram. For 
the AGN where this correction failed (i.e. L[ iii],agn > ^ioiii],obs), we set 1^[oiii],agn = 
L[oni].obs- For the composite and star-forming galaxies where this correction failed (i.e. 
^[oiii],agn < 0), we set L[ iii],agn = 0. A vast majority of the star-forming galaxies had 
L[oiii],agn set to 0. 

Using the sources that had [Nell] detections, we calculated a and (3 by using a least 
squares trimmed (LTS) regression, a robust regression technique that removes outliers, on 
the above equation. As each execution of the LTS regression results in different values for 
a and (3, we performed this calculation 1000 times to determine the average values and 
spreads on the distributions. The resulting values, using a 3a clipping mean, are a=89±l 
and /9=111±7. We note that the difference between these values and using a straight average 
are negligible for a and small for (c.f. 109±15). In Figure PT4"l we plot Lmir vs. oLsfr + 
(3Lagn (i-e. 89 x L[ NeII j + 124 x L[ iii],agn), with the overplotted line indicating where 
the two quantities are equal. No systematic offsets among galaxy sub-type are present. This 
relationship can therefore be valuable for estimating the MIR flux due to star- format ion (or 
conversely from the torus) in active galaxies. 

To test this decomposition, we have plotted axL[ NeII }/LMiR vs. the PAH EW at 7.7/iin 
in Figure [151 As we have demonstrated above, this PAH macro- molecule/small-grain is the 
least affected by AGN activity (i.e. not suppressed like the 11. 3/im feature nor mildly boosted 
like the 17/im feature). Hence, to first order the EW of the 7.7/im PAH feature is tracing the 
"extra" MIR emission from AGN heated dust. If our decomposition is correct, the relative 
amount of MIR emission due to star formation activity, i.e. axL^ e jj]/LMiR, should increase 
monotonically with increasing EW of the PAH 7.7/im feature. The dotted dashed lines have 
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unity slope and were chosen to intersect the points where the MIR emission is solely due 
to star formation (e.g. log (ax L[ NeII j/L M j R ) = dex) and fiducial PAH EW values that 
approximately enclose the star-forming galaxy locus. Many of the AGN fall within the region 
demarcated by these two lines. The dashed line is the best-fit line to the AGN sources where 
the slope is held constant at unity. This fit neatly bisects the above region. These results 
suggest that our empirical decomposition works reasonably well in describing the relative 
contributions of star formation and AGN activity to the MIR emission. 



7. Conclusions 

We have explored the infrared and optical parameter space where star-forming galaxies, 
composites and AGN live to analyze diagnostics that parameterize host galaxy star for- 
mation. Using star-forming galaxies as a control sample, we have investigated which star 
formation diagnostics are least affected by AGN activity. We have studied diagnostics that 
trace the interplay between AGN and starburst activity and tested whether a smooth transi- 
tion exists over a range of radiation field hardness. Finally, using the results of this analysis, 
we present an empirical decomposition of the MIR flux into a star-forming and an AGN 
component. Our overall result is that the optical and mid-IR diagnostics of star formation 
and of the relative importance of young stars and the AGN generally agree well. Our more 
specific results are summarized as follows: 

SFR diagnostics The SDSS derived SFR/jf, er and the [Nell] 12.81/xm luminosity agree 
well, and are the most reliable SF proxies we have considered for use in star-forming 
galaxies, composites and AGN. The sum of [Nell] and [Nelll] 15.56/im is systemati- 
cally offset to higher values for AGN due to the active nucleus boosting the [Nelll] 
flux. The aggregate PAH luminosity {LT,PAH7.7,i\.z,nnm) relative to L[jvejj] is sup- 
pressed in AGN compared with both quiescent and composite galaxies. Though the 
(LsPAJ3"7.7,n.3,i7/im)/SFRjj6 er . ratio agrees between star-forming galaxies and Sy2s, anal- 
ysis of PAH flux ratios indicate abnormal behavior in some AGN dominated systems, 
such as suppression of the 11.3/im feature, suggesting diminished PAH emission in 
strong AGN rather than [Nell] enhancement. The disparate slopes of (LspAH7.7,n.3,i7>m) /SFR/j& er 
between Sy2s and star-forming galaxies also hint that AGN may contaminate the emis- 
sion from PAH macro-molecules. 



Ionization Field Hardness We have used the optical D parameter, which indicates 
the distance a source l ies from the star-fo rming galaxy locus on the BPT diagram 
( iKauffmann et al.ll2003l ). and L[jv e y]/L[Are/7] ( jArmus et al.ll2007l ; |Pereira-Santaella et al. 
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20101 ) as probes of the incident radiation field hardness. The infrared diagnostics that 
parameterize the relative contributions of AGN to star-forming activity that we have 
considered (tt2o-30/im ; PAH EWs and PAH flux ratios) only show significant trends 
with ionization field hardness in the AGN population. This result suggests that rather 
than parameterizing incident radiation hardness in general, these IR diagnostics are 
only effective descriptors of ionization field hardness in AGN. 

a 2 o-30/im In AGN, the spectral slope between 20-30 /xm (measured by the spectral 
index o^n-a rym) has been shown to steepen in the p resence of cold dust due to star 



formation ([Buchanan et al.ll2006l ; iDeo et al.l 120091 ). For the AGN in our sample, 
«2o-30/im is significantly anti-correlated with proxies of ionization field hardness, 
albeit with wide scatter, however this result does not hold when expanded into 
the regime of quiescent star forming galaxies. 

PAH EW PAH EWs are significantly anti-correlated with the increasing hardness of 
radiation field, but this is only true when considering the AGN. When considering 
just the star-forming galaxies, no trend exists between incident radiation hardness 



(2009 


) and 


Trever et al. 


(2010) 



hypothesis that the mid-IR continuum has a significant contribution from dust 
heated by the AGN only in objects where the ionizing radiation field is dominated 
by the AGN. For the 11.3/xm PAH feature, our analysis of the PAH flux ratios 
suggests that destruction of the PAH macro-molecule in strong AGN could also 
contribute to suppression of the PAH 11.3/xm EW. 

PAH flux ratios We find that the luminosity of PAH 11.3 \xm feature is signifcantly 
suppressed relative to the luminosities of the PAH 7.7 fim and 17 \xm features 
and relative to the optical and mid-IR derived star formation rates in the most 
AGN-dominated systems. Thus, some care is required in using the PAH 11.3/iin 
luminosity as a proxy for the SFR in AGN-dominated systems. 

Composite Systems Galaxies that are optically classified as composite systems are more 
akin to quiescent star forming galaxies in terms of the mid-IR parameters that trace 
star formation: the smaller relative energetic importance of the AGN in composites 
does not seem to affect [Nelll] and PAH luminosities. 



[NeV] as unambiguous AGN signature Due to the high ionization p otential of 



(97.1 eV), it is poten tially an unequivocal signature of AGN activity ( IGilli et al. 



Goulding et al. 



NeV] 



2010 



20101 ) . One galaxy that is optically classied as a quiescent star-forming 



system in our sample has a [NeV] detection, yet other IR parameters (i.e. PAH EW 
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values) are normal for a pure star forming galaxy. The nature of this galaxy as a strong 
"hidden" AGN is not clear. 



Disentangling MIR Emission Using multiple linear regression, we fit the relation Lmir 
= aLsFR + (3I->agn to our full sample to determine a and (3. We set Lsfr = L[ NeII } and 
^agn=^\oiii],agn where in L[ iii],agn, we h ave subtract e d out the estimated starburst 
contribution to the [OIII] flux using Eq. 3 in Iwild et all fcoioh . We find a=89±l and 
/3=111±7. This decomposition can be useful in estimating the MIR emission due to 
the circumnuclear AGN obscuration or conversely host galaxy star formation in AGN. 
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Fig. 1. — Left: Ljjve/j] vs - SDSS SFR/jb er . Right: Distribution of log(L[7v e //]/SDSS SFR/ife er ). 
The symbol coding is as follows: S5 - crosses, SSGSS - circles, [OIII] sample - squares, 
12/im sample - diamonds. Blue data points represent star-forming galaxies, green indicates 
composites and red denotes AGN. The dashed-dotted line indicates the best fit to the star- 
forming galaxy sub-sample and the red dashed line shows the best fit to the AGN population. 
Active galaxies agree with quiescent in terms of the L^ Nen ^ distribution relative to SDSS 

SFR fib er . 
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Fig. 2. — Left: L[jve/r] + L[mj/] vs - SDSS SFRj iber . Right: Distribution of 
log((L[jv e //]+L[Ar e /7/]/SDSS SFR/jfe er ). We note that the relatively large number of non- 
detections causes the best-fit relation from survival analysis to lie below the majority of 
plotted upper limits. Nevertheless, the AGN clearly have boosted L[ NeII } + h^ NeIII ] values 
compared with star-forming and composite galaxies, which can be seen by comparing the 
best-fit trend to the Sy2s (red dashed line) with the star-forming galaxies (blue dotted- 
dashed line). Compared with Figure 1, a clear offset is discernable between the best-fit 
relations between Sy2s and quiescent galaxies, suggesting AGN enhances [Nelll] emission. 
Color, symbol and line coding same as Figure 1. 
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Fig. 3. — Top Left: LEPAH7.7,ii.3,i7/im v s SDSS SFR/jft er . Top Right: Distribution of 
log(LEPAi?7.7,ii.3,i7/um/SDSS SFRj iber ). Bottom Left: 1^t,pah7. 7, 11.^,17 vs L[7v e //]. Bottom 
Right: Distribution oi\og{\jY,pAH7.7,ii.^i7^m/^[Neii])- The sum of the PAH luminosities rela- 
tive the SFR/jher is consistent among star-forming galaxies, composites and AGN. However, 
when compared with L[jveJj], the PAH emission in AGN is suppressed by about 32% on- 
average in relation to the observed PAH emission in star-forming galaxies and composites: 
the mean log(Lsp J 4^7.7,ii.3,i7/im)/L[Ar e //] is 1.50 dex for star-forming galaxies, yet 1.33 dex 
for Sy2s. The best-fit trend to the Sy2s (red dashed line) also illustrates the different re- 
lationship between the PAH luminosities in AGN compared with quiescent galaxies (blue 
dotted-dashed line). Color, symbol and line coding same as Figure 1. 
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Fig . 4. — Left: BPT dia gram for all s ources in our sample. Th e unbroken black lines indicate 
the lKewley et al.l ( l200ll ) (upper) and iKauffmann et al.l (120031 ) (lower) demarcations between 
star-forming galaxies, composites and AGN. The dashed black lines illustrate the dividing 
lines we used to separate the sample to stack the short-high IRS spectra. Symbol and color 
coding is same as Figure 1. Right: After removing the low-metallicity systems and using the 
demarcations shown in the left-hand plot, the groups into which our sample were separated 
our indicated by the colors in this plot: star-forming group 1 - dark blue, star-forming group 
2 - light blue, composite group 1 - dark green, composite group 2 - light green, AGN group 
1 - gold, AGN group 2 - red. 
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Fig. 5. — Stacked SH spectra around the [NeV] emission line at 14.32 /jm for each group of 
galaxies illustrated in Figure H]b). The dashed line indicates the nominal [NeV] wavelength. 
[NeV] was detected above the 5a level in AGN group 2 only. 
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Fig. 6. — Left: L[ NeV }/L[ NeII ] vs. the optical D parameter. Right: \j\ i N e v}/^-'[Neii} vs- 
a 2o-30^tm- These quantities are statistically correlated and anti-correlated: p=0.562 and 
p=-0.568, respectively. Color and symbol coding same as Figure SJ The red triangle indi- 
cates the [NeV] detection for the stacked spectra in AGN group 2. For these and subsequent 
plots with IR emission line and PAH fluxes, error bars are smaller or on the order of the 
symbol size and therefore not plotted. 
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Fig. 7. — a) L[Ar e y]/L[Ar e //] vs. PAH EW at 7.7 fim. b) L[7v e y]/L[7v e i7] vs. PAH EW at 
11.3 fim. c) L[jveV]/L[jVe//] vs - PAH EW at 17 /im. These quantities are statistically anti- 
correlated: p=-0.702, p=-0.779 and p=-0.628, respectively. The red triangle indicates the 
[NeV] detection for the stacked spectra in AGN group 2. Color and symbol coding same as 
Figure 1. 
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Fig. 8. — a2o-30/«n vs. the optical D parameter. These parameters are correlated for AGN, 
p=-0.474, but not for the full sample, p=-0.195. Color and symbol coding same as Figure HJ 
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Fig. 9.— The optical D parameter vs a) PAH EW at 7.7 fim, b) PAH EW at 11.3 pm and 
c) PAH EW at 17 pm. These quantities are statistically anti-correlated: p=-0.679, p=-0.747 
and p=-0.658, respectively. Color and symbol coding same as Figure 1. 
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Fig. 10.— L [NeV] /L [NeII] vs. a) L PAH7 , 7fim /L PAH17fim (no significant trend, p=-0.241), b) 
Lp4ii-ii.3/xm/LpAifi7/xm (significant anti-correlation, p=-0.564) and c) L PAH7 , 7fMm /L PAH11 , 3fim 
(significant correlation, p=0.611). The red triangle indicates the [NeV] detection for the 
stacked spectra in AGN group 2. Color and symbol coding same as Figure 1. 
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Fig. 11— Optical D parameter vs. a) L PAH7 , 7fim /L PAH17fim , b) L PAH11 , 3fim /L PAH17fim , p=- 
0.364) and c) l^pAH7.7 i im/^- , PAHix.z IJ ,m- Color and symbol coding same as Figure 1. 
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Fig. 12. — The optical D parameter vs. PAH EWs. Imposing a redshift cut, to test for 
possible systematic differences due to physical sizes probed between closeby and far away 
sources, does not alter the trends seen in the full sample. Color and symbol coding same as 
Figure 1. 
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Fig. 13. — For sources at z> 0.03, optical D parameter vs. PAH ratios shown in Figure [TBI 
Imposing a redshift cut does not alter the trends seen in the full sample. Color and symbol 
coding same as Figure 1. 
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Fig. 14 — L M ir vs. aLsFR + P^agn, where L S fr = L[Neii] and L AG n = L[oiii],agn- a and 
(3 were found by using a least trimmed squares regression: a=89±l, /3=111±7. The dashed 
line indicates where the two quantities are equal. Color coding same as Figure 1. 
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Fig. 15. — Log (axL[ NeII }/LMiR) vs. Log (PAH EW) at 7.7/xm. The two dotted-dashed lines 
have unity slope and intercept points where the MIR emission is due purely to star formation 
according to our decomposition (i.e. ax L[jv e jj] = Lmir) at two fiducial values of PAH 
7.7/im EW (1.0 dex and 0.5 dex) chosen to bracket the star- forming galaxies. The dashed 
line indicates the best-fit through the AGN, with the slope held constant at unity. The AGN 
generally fall within the region indicated by the lines, which is consistent with a monotonic 
decrease in PAH EW, denoting the "extra" MIR emission due to AGN, correlating with 
the empirial MIR AGN parameterization (i.e. decreasing ax\j\ iNeII }/~L MIR values indicate 
greater AGN contribution to the MIR luminosity). Color and symbol coding same as Figure 
1. 



Table 1. Emission Line Luminosities for SSGSS Sample 



Source Number 


RA 


Dec 


z 


Type 


[Nell] 12.81/im 


[NeV] 14.32//m 


[Nelll] 15.56/xm 












10 41 erg/s 


10 41 erg/s 


10 41 erg/s 


33 


160.291 


56.932 


0.185 


AGN 


8.12± 0.70 




3.34± 0.57 


61 


160.952 


58.197 


0.073 


AGN 


3.19± 0.08 




0.30± 0.05 


11 


161.720 


56.252 


0.047 


Comp 


0.17± 0.01 




< 0.06 


12 


162.268 


56.224 


0.072 


Comp 


3.04± 0.09 




0.31± 0.04 


15 


161.758 


56.307 


0.153 


Comp 


4.40± 0.30 




< 1.64 



1 This table is available in its entirety in machine readable format on-line. We list several row here as a 
guide. 



Table 2. Emission Line Luminosities for S5 Sample 



RA 


Dec 


z 


Type 


[Nell] 12.81/xm 


[NeV] 14.32//m 


[Nelll] 15.56/im 










10 41 erg/s 


10 41 erg/s 


10 41 erg/s 


36.380 


-8.419 


0.055 


AGN 


1.16± 0.05 




1.02± 0.07 


36.435 


-7.868 


0.077 


AGN 


2.72± 0.12 




1.03± 0.08 


132.186 


53.746 


0.083 


AGN 


4.72± 0.14 


1.18± 0.09 


2.54± 0.08 


132.714 


54.823 


0.081 


AGN 


3.17± 0.11 




2.25± 0.07 


132.815 


54.986 


0.083 


AGN 


3.64± 0.16 




1.68± 0.09 



1 This table is available in its entirety in machine readable format on-line. We list several 
row here as a guide. 



Table 3. Emission Line Luminosities for [OIII] Sample 



Galaxy 




RA 


Dec 


z 


[Nell] 12.81/im 


[NeV] 14.32 / um 


[Nelll] 15.56/im 












10 41 erg/s 


10 41 erg/s 


10 41 erg/s 


NGC 291 




13.375 


-8.768 


0.019 


1.61± 0.04 


0.60± 0.02 


1.12± 0.02 


Mrk 609 




51.356 


-6.144 


0.034 


5.51± 0.12 




1.80± 0.06 


IC 0486 




120.087 


26.614 


0.027 


1.10± 0.10 


0.55± 0.06 


1.04± 0.05 


2MASX J08035923- 


h2345201 


120.997 


23.756 


0.029 




< 0.14 


< 0.42 


2MASX J08244333- 


F2959238 


126.180 


29.990 


0.025 


0.54± 0.06 


0.55± 0.07 


0.70± 0.06 



1 This table is available in its entirety in machine readable format on-line. We list several row here as a guide. 



Table 4. Emission Line Luminosities for 12/xm Sample 



Galaxy 


RA 


Dec 


z 


Type 


[Nell] 12.81//m 


[NeV] 14.32/im 


[Nelll] 15.56/im 












10 40 erg/s 


10 40 erg/s 


10 40 erg/s 


IRAS 00198-7926 


5.473 


-79.169 


0.073 


AGN 


35.65± 1.45 


49.05± 1.53 


60.06± 1.31 


NGC 0424 


17.865 


-38.083 


0.012 


AGN 


3.06± 0.19 


4.64± 0.27 


5.67± 0.17 


NGC 1068 


40.670 


-0.013 


0.004 


AGN 


5.34± 0.28 


10.11± 0.48 


15.02± 0.32 


NGC 1144 


43.790 


-0.178 


0.028 


AGN 


10.79± 0.26 


0.61± 0.10 


2.88± 0.05 


NGC 1320 


51.203 


-3.042 


0.009 


AGN 


0.47± 0.02 


0.59± 0.02 


0.78± 0.02 



1 This table is available in its entirety in machine readable format on-line. We list several row here as a 
guide. 



Table 5. PAH Luminosities and EWs for SSGSS Sample 



Src Num 


Type 


L7.7/itm 


EW 7.7/mi 


Lll.3/im 


EW 11.3/xm 


L 17 fim 


EW 17/im 






10 41 erg/s 


/im 


10 41 erg/s 


/im 


10 41 erg/s 


/im 


13 


AGN 


16.59± 0.82 


1.23 


7.61± 0.32 


0.39 


6.12± 1.09 


0.27 


21 


AGN 


30.66± 1.73 


2.95 


10.80± 0.60 


1.43 


6.44± 0.49 


1.09 


23 


AGN 














33 


AGN 


157.18± 6.51 


3.14 


53.40± 1.08 


1.01 


37.98± 3.62 


0.66 



50 AGN 



1 This table is available in its entirety in machine readable format on-line. We list several row here as a 
guide. 



Table 6. PAH Luminosities and EWs for S5 Sample 



RA 


Dec 


Type 




EW 7.7/im 


Lll.3£tm 


EW 11.3/xm 




EW 17/xm 








10 41 erg/s 


/im 


10 41 erg/s 


fxm 


10 41 erg/s 


/im 


36.380 


-8.419 


AGN 


24.01± 2.34 


5.74 


7.92± 0.71 


3.62 


15.13± 2.76 


2.28 


36.435 


-7.868 


AGN 


59.59± 4.64 


1.44 


16.41± 0.62 


0.53 


7.79± 1.54 


0.28 


132.186 


53.746 


AGN 


144.14± 5.69 


4.81 


41.27± 0.81 


1.81 


32.40± 1.82 


1.38 


132.714 


54.823 


AGN 


54.59± 5.05 


4.51 


21.61± 0.71 


1.89 


16.96± 1.69 


1.04 


132.815 


54.986 


AGN 


107.85± 3.30 


5.05 


29.54± 0.49 


1.62 


19.68± 1.23 


1.17 



This table is available in its entirety in machine readable format on-line. We list several row here as a guide. 



Table 7. PAH Luminosities and EWs for [OIII] Sample 



Galaxy 




EW 7.7 nm 


Lll.3^m 


EW 11.3/im 




EW 17/xm 




10 41 erg/s 


/im 


10 41 erg/s 


/im 


10 41 erg/s 


/im 


NGC 0291 


34.69± 0.95 


10.07 


7.44± 0.15 


0.75 


9.20± 0.15 


0.50 


Mrk 0609 






52.35± 0.78 


1.08 


54.58± 0.54 


0.94 


IC 0486 


22.69± 0.47 


0.75 


4.15± 0.19 


0.12 


5.77± 0.12 


0.16 


2MASX J08035923+2345201 






1.49± 0.05 


0.30 


2.10± 0.13 


0.40 


2MASX J08244333+2959238 






2.80± 0.16 


0.04 


5.67± 0.29 


0.10 



1 This table is available in its entirety in machine readable format on-line. We list several row here as a guide. 



Table 8. PAH Luminosities and EWs for 12 /jm Sample 



Galaxy 


Type 




EW 7. 7 pan 


Lll.3Atm 


EW 11.3/xm 


Ll7 /j,m 


EW 17 pan 






10 41 erg/s 


pan 


10 41 erg/s 


/im 


10 41 erg/s 


/im 


NGC 0424 


AGN 


10.25± 0.65 


0.09 


4.46± 0.35 


0.06 


8.05± 0.61 


0.19 


NGC 1144 


AGN 


258.85± 1.37 


2.21 


92.63± 0.99 


1.83 


51.43± 0.73 


1.56 


NGC 1320 


AGN 


9.69± 0.38 


0.52 


2.45± 0.08 


0.16 


1.31± 0.12 


0.11 


NGC 1386 


AGN 


0.84± 0.02 


0.71 


0.16± 0.00 


0.14 


0.26± 0.01 


0.27 


NGC 1667 


AGN 


61.97± 1.14 


3.85 


18.76± 0.22 


2.18 


9.18± 0.13 


1.63 



1 This table is available in its entirety in machine readable format on-line. We list several row here as a 
guide. 



Table 9. Ratios 1 of SFR proxies 



Ratio SF 2 N 2 Composites N AGN N Comments 



L[AT e 77]/SFR/jb er 


40.97j£| 


223 3 


40 Q1+ - 08 


32 


40.92lg|! 


39 3 




(L [Nell] +L [Nelll] ) / SFRjj feer 


40.841°;^ 


223 3 


4fl 74+ - 05 


32 


41 11+ - 21 


39 




LsPAH7.7,11.3,17Atm/SFRjjfe er 




261 3 




46 


49 4Q+ - 16 


35 3 




LEPAH7.7,11.3,17Atm/L[7V e //] 


1.50 4 


207 3 


i.oy_ 22 


31 


1 oq+0.15 
i -°°-0.59 


21 


PAH ULs 


L[Ar e 7/]/LEPAi?7.7,11.3,17Atm 


-1.61 4 


141 3 


-1 fi2 +a04 


25 


-1 50+ 002 

l.OU_ 21 


17 


[Nell] ULs 



1 Ratios in log space. Ratios were calculated using survival analysis. Error bars represent 
the deviation from the values at the 75th and 25th percentiles from the mean. 

2 "SF" refers to star-forming galaxies. 

3 Number of sources in each galaxy sub-group used in calculation of mean values. 

4 Upper limit data point changed to a detection for Kaplan-Meier computation, which biases 
the mean estimate. 

4 Error in percentile calculation, so error bars are not reported. 



Table 10. Two-Sample Tests on SFR Proxies between Galaxy Sub-Samples 1 



Ratio 


SF vs. Composites 
P 


N 2 


Composites vs. AGN 
P 


N 


SF vs. AGN 
P 


N 


Comments 3 


L[7Ve77]/SFRj-j ber 


0.127 


255 


0.598 


71 


0.624 


223 




(L[JVeJJ] +L [Nell I] ) / SFRj i6er . 


0.175 


255 


0.0005 


71 


<lxl0~ 4 


262 




LsPAH7.7,11.3,17Atm/SFR/jb er . 


0.331 


307 


0.365 


81 


0.758 


296 




LsPAff7.7,11.3,17/im/L[jVelJ] 


0.213 


238 


0.0002 


77 


0.003 


253 


PAH ULs 


L[JVe//]/L£PAH7.7,11.3,17/xm 


0.945 


166 


0.006 


42 


0.0004 


158 


[Nell] ULs 



1 "SF" refers to star-forming galaxies. 
2 Number of sources in two-sample tests. 

3 We report the survival analysis probability from Gehan's Generalized Wilcoxon Test, permutation variance. For 
the PAH and [Nell] ratios, survival analysis was performed twice: once to account for [Nell] upper limits and once to 
accommodate the PAH luminosity upper limits. A P value <0.05 indicates that the two sub-populations differ at a 
statistically significant level. 
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Table 11. Correlation between Diagnostics of Relative Importance of AGN Activity to 

Star formation 1 



Diagnostic 


P 2 


P 3 

1 uncorr 


L[jVeV]/L[iVe/7] V S- D 


0.562 


4xl0" 4 


L[Ar e y]/L[Ar e //] VS. C*20-30/jm 


-0.568 


6xl0" 4 


L [NeV] /L [NeII] vs. PAH EW 7.7 //in 


-0.702 


lxl0~ 4 


L[JVeV]/L[ATe77] VS. PAH EW 11.3 /im 


-0.779 


< 1 x 10~ 4 


L [iVe y]/L[jve//] vs. PAH EW 17 jim 


-0.628 


2xl0" 4 



1 Used ASURV to account for upper limits in L[Ar e y]- 

2 p is the correlation coefficient from which ~P uncoTr is de- 
rived. 

Probability the two quantities are not correlated. A 
Puncorr value <0.05 indicates a significant correlation. 
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Table 12. Relationship between PAH EWs and Ionization Field Hardness 



Diagnostic 


P 1 


N 2 


D vs. PAH EW 7.7 /im 


-0.679 


344 


D vs. PAH EW 11.3 /im 


-0.747 


368 


D vs. PAH EW 17 /im 


-0.658 


281 


Star-forming Galaxies and Composites Sub-sample 


D vs. PAH EW 7.7 /im 


-0.227 


291 


D vs. PAH EW 11.3 /im 


-0.125 


303 


D vs. PAH EW 17 /im 


-0.230 


221 


Sources with z > 0.03 


D vs. PAH EW 7.7 /im 


-0.566 


318 


D vs. PAH EW 11.3 /im 


-0.653 


336 


D vs. PAH EW 17 /im 


-0.549 


249 



l p is the correlation coefficient. Higher \p\ values 
indicate more significant (anti-)correlations. 

2 Number of sources. 
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Table 13. Relationship between PAH Ratios and Ionization Field Hardness 



Diagnostic 


P 1 


N 2 


[NeV]/[NeII] 3 vs.: 


L p AH7. /L p AHl 


-0.241 


33 


LpAffll.3/im/LpAH17/im 


-0.564 


37 


LpAP"7.7Mm/LpAP11.3/im 


0.611 


34 


D 4 vs.: 


LpAH7.7/xm/LpA/n7/xm 


-0.193 


266 


LpAH11.3nm/I J PAH17iJ,m 


-0.364 


278 


LpAff7.7/xm/LpA.ffll.3jim 


0.179 


340 


AGN sub-sample, D 4 vs: 


LpAP7.7Mm/LpAP17Mm 


0.082 


47 


LpAP11.3Atm/LpAH17Aim 


-0.247 


58 


LpAP7.7^m/LpAP11.3Atm 


0.431 


51 


Sources with z > 0.03, 


D 4 vs.: 




LpAP7.7Mm/LpAP17Mm 


-0.283 


240 


LpAi?11.3/im/LpAH17/im 


-0.356 


246 


^PAH7.7fMm/^PAH11.3fim 


0.084 


314 


AGN with z> 0.03, D 4 


vs: 




LpAP7.7^m/LpAP17^m 


-0.208 


22 


Lp Aif" 1 1 . 3/«n / L p AH1 7/im 


-0.554 


27 


LpAf/7.7^m/LpAP11.3Atm 


0.590 


26 
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1 p is the correlation coefficient. 
Higher \p\ values indicate more signifi- 
cant (anti-) correlations. 

2 Number of sources. 

3 p calculated from survival analysis. 

4 p calculated from linear regression. 



